Segregation in hard sphere mixtures under gravity. An extension of Edwards 
approach with two thermodynamical parameters 
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We study segregation patterns in a hard sphere binary 
model under gravity subject to sequences of taps. We discuss 
the appearance of the "Brazil nut" effect (where large parti- 
cles move up) and the "reverse Brazil nut" effects in the sta- 
tionary states reached by "tap" dynamics. In particular, we 
show that the stationary state depends only on two thermody- 
namical quantities: the gravitational energy of the first and of 
the second species and not on the sample history. To describe 
the properties of the system, we generalize Edwards' approach 
by introducing a canonical distribution characterized by two 
configurational temperatures, conjugate to the energies of the 
two species. This is supported by Monte Carlo calculations 
showing that the average of several quantities over the tap dy- 
namics and over such distribution coincide. The segregation 
problem can then be understood as an equilibrium statistical 
mechanics problem with two control parameters. 

Segregation is an ubiquitous and intriguing phe- 
nomenon observed in vibrated granular mixtures such as 
powders or sand: in presence of shaking the system is not 
randomized, but its components tend to separate [jl],^). 
An example is the so called "Brazil nut" effect (BNE) 
where, under shaking, large particles rise to the top and 
small particles move to the bottom of the container. Sev- 
eral mechanisms have been proposed to explain these 
phenomena, as for instance "percolation" H, geomet- 
rical reorganization , convection Q , "condensation" 
or inertia ||. Interestingly, by changing grains sizes 
or mass ratio a crossover towards a "reverse Brazil nut" 
effect (RBNE) was more recently discovered |?J where 
small particles segregates to the top and large particles 
to the bottom. On the whole the criteria to predict segre- 
gation in a mixture, an issue of large practical relevance, 
are still largely unknown |jj . 

In the perspective of the Statistical Mechanics ap- 
proach originally proposed by Edwards to describe pow- 
ders ||, we study here segregation patterns, in absence 
of convection mechanisms, of a binary hard sphere model 
under gravity subject to sequences of "taps" . We show 
that a Statistical Mechanics description of segregation is 
indeed possible by the introduction of two "configura- 
tional temperatures" which characterize the macroscopic 
status of the system. 

In Edwards' approach granular media at rest (i.e., 
when grains have no kinetic energy and, thus, the bath 
temperature is zero) can be described in terms of Statis- 
tical Mechanics where a single control parameter, called 
compactivity, plays the role of the temperature of usual 
thermal systems pIMpTlfflO . We show that the sta- 



tionary states of the present mixture are indeed indepen- 
dent on the sample history as in a "thermodynamics" 
system, but are characterized by two control parameters, 
such as the gravitational energy of the first and of the 
second species. We postulate that the probability dis- 
tribution of the system in the stationary state can be 
obtained by the principle of maximum entropy under the 
constraint that the two energies are independently fixed. 
This leads to a Gibbs canonical distribution character- 
ized by two configurational temperatures, conjugate to 
the energies of the two species, generalizing in this way 
Edwards approach [llj . The postulate is supported by 
Monte Carlo simulations showing for several quantities 
that statistical averages over the tap dynamics coincide 
with those over such an ensemble distribution. The seg- 
regation problem can then be studied as standard Statis- 
tical Mechanics problems using such a distribution and, 
for instance, a segregation "phase" diagram can be the- 
oretically derived. 

We consider a binary hard sphere system made of two 
species 1 (small) and 2 (large) with grain diameters ao 
and V2ao, under gravity on a cubic lattice of spacing 
ao = 1 (see inset of FigJ|) . The hard core potential Tihc 
is such that two nearest neighbor particles cannot over- 
lap. This implies that only pairs of small particles can 
be nearest neighbors on the lattice. The system overall 
Hamiltonian is: 



H = TLhc + niigHx + m 2 gH 2 



(1) 



where H\ = J^j z i an d H 2 = Y^T' Zi i-> ^ ne height of site 
i is Zi and the two sums are over all particles of species 
I and 2 respectively. We set the units such that the two 
kinds of grain have masses mi — I and m 2 — 2, and the 
gravity acceleration is g = 1. In the above units, the 
gravitational energies in a given configuration are thus 
E 1 = Hi and E 2 = 2H 2 . 

Grains are confined in a box of linear size L between 
hard walls and periodic boundary conditions in the hori- 
zontal directions. The system is subject to a Monte Carlo 
dynamics made of sequences of "taps" ]T^ , ^0| , pTt . Dur- 
ing each "tap" the system evolves for a time tq (the "tap 
duration" , in units of lattice sweeps) in presence of a 
finite thermal bath of temperature TrS afterwards it is 
suddenly frozen at zero temperature in one of its stable 
states. These states correspond to a local minimum or 
saddle of the energy, such that any particle movement 
does not decrease the energy, and have been called "in- 
herent states" in analogy with the terminology of 
glasses IPTy . After each tap, when the system is at rest, 



,(2) 
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we record the quantity of interest, which are a function 
of the number of "taps" , t. 

Ni grains of type 1 and N2 grains of type 2 are ini- 
tially prepared in a random loose stable pack (l8) . Under 
"shaking" the system is let to approach a stationary state 
for each value of the tap parameters Tr and To, as shown 
in Fig-0 for the large grains height /12 (t) — Hi (t) /N2 ■ In 
Fig.g, we plot as function of Tr (for several values of To) 
the asymptotic value of the vertical segregation parame- 
ter, i.e., the difference of the average heights of the small 
and large grains at stationarity, Ah(Tp,r ) = h± — h 2 - 
Here hi and hi are the average of Hi/N x and H2/N2 
over the tap dynamics in the stationary state. 

Fig.|| shows that at high Tr the Brazil Nut Effect 
(BNE) of granular media is observed: "large" grains are 
found on the top of the system (i.e., Ah < 0). For 
small values of Tr, below a crossover amplitude, Tp(ro), 
a smooth crossover to the Reverse Brazil Nut Effect 
(RBNE) is instead recorded, with Ah > 0. The RBNE 
was hrst discovered in recent Molecular Dynamics simu- 
lations . Our aim below is to show that the properties 
of these segregation patterns can be explained in terms 
of a statistical mechanics approach to the problem. 

The results given in the main panel of Fig. |^ apparently 
show that Tr is not a right "thermodynamical" parame- 
ter, since sequences of "taps" with different r give dif- 
ferent values for the system observables. However, if the 
stationary states corresponding to different tap dynam- 
ics (i.e., different Tr and To), are indeed characterized 
by a single "thermodynamical" parameter a la Edwards, 
the curves of Fig.|| should collapse on an universal mas- 
ter function when Ah(Tr, tq) is parametrically plotted as 
function of an other macroscopic observable such as the 
average energy, e(T r ,r ) = (E 1 +E 2 )/N = ^hi + ^h 2 
(where N = N\ + N 2 ). This collapse of data, actually 
found in other situations JlO| , pT[ , is clearly not observed 
here, as it is apparent in the inset of Fig.^J (see also 

We show, instead, that two macroscopic quantities, 
such as hi and /12, may be sufficient to characterize 
uniquely the stationary state of the system. Namely, we 
show that any macroscopic quantity Q, averaged over 
the tap dynamics in the stationary state, is only de- 
pendent on hi and /12, i.e., Q — Q(hi,h2). We have 
checked that this is the case for several independent ob- 
servables, such as the number of contacts between large 
particles, N c , the density of small and large particles on 
the bottom layer, p\ and p\, and others. As shown in 
Fig.|3|, we find, in particular, with good approximation 
that: iV c ~ N c (e) = N c (^hi + ^h 2 ), p\ * ^(fc 2 ), 
p\ ~ p\(hi). Therefore we need both hi and h% (or 
equivalently Ei and E2) to characterize unambiguously 
the state of the system independently on the previous 
history (i.e., in our case independently on the particular 
tapping parameters Tr and To). 

In this respect the stationary state can be genuinely 
considered as a "thermodynamical state", and conse- 
quently one can attempt to construct an equilibrium sta- 



tistical mechanics, as originally suggested by Edwards. 
This will allow to calculate any quantity in the station- 
ary state using an appropriate ensemble, irrespective of 
the particular history. 

In the present paper we generalize Edwards' approach 
by introducing a canonical distribution characterized by 
two configurational temperatures, conjugate to the ener- 
gies of the two species, and verify by Monte Carlo calcu- 
lations that the average of several quantities over the tap 
dynamics and over such distribution coincide. More pre- 
cisely we have to find the probability distribution, P r , of 
observing the system in the generic inherent state r cor- 
responding to an energy E Xr for the small particles and 
E2r for the large particles (see To answer this, we 

assume that the microscopic distribution is given by the 
principle of maximum entropy S — — ^ r P r In P r under 
the condition that the average energy Ei = ^ r P r Ei r 
and E2 = Y^ ir P r E2r are independently fixed. This can 
be done by introducing two Lagrange multipliers /3i and 
p2i which are determined by the constraint on Ei and 
E2, and can be considered as the "inverse configurational 
temperature" of species 1 and 2 |nj. This procedure 
leads to the Gibbs result: 

Pr = ~ ^ (2) 

where Z = ^ r exp [— j3\E\ r — /3a-E2r] and, in the ther- 
modynamic limit, the entropy, S, is given by: 

S = ]nSl(E 1) E 2 ), (3) 

and Pi and 02: 

dlnn(Ei,E 2 ) d\nn{Ei,E 2 ) ... 

Pl = — dEi — ' 02 = — m — ■ (4) 

Here £l(Ei,E 2 ) is the number of inherent states corre- 
sponding to energy E\ and i?2- 

The hypothesis that the ensemble distribution at sta- 
tionarity is given by eq.(||) can be tested as follows. We 
have to check that the time average of any quantity, 
Q{hi,%2), as recorded during the taps sequences in a 
stationary state characterized by given values hi and h 2 , 
must coincide with the ensemble average, (Q)(hi, /12), 
over the distribution eq.(|^). To this aim, we have calcu- 
lated the ensemble averages (N c ), (p\), (p\), for different 
values of 0i and P2] then we have expressed parametri- 
cally (iV c ), (py), (p\), as function of the average of hi and 
/i2, and compared them with the corresponding quanti- 
ties, N c , p\ and p\, averaged over the tap dynamics. The 
two sets of data are plotted in FigJ| showing a rather 
good agreement (notice, there are no adjustable param- 
eters). Technically, in order to calculate the ensemble 
averages we simulated, with standard Monte Carlo meth- 
ods, the model with Tt from eq.(Q) where we impose that 
the only accessible states are the inherent states. This 
can be done by using a configurational Hamiltonian 

Hhc + PxmigHi + /3 2 m2gH 2 - Inn, (5) 
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where II — 1 or depending whether the configuration is 
an inherent state or not. In this way the weight e _Wco,l - f 
gives the ensamble distribution (||) (for further details see 
jpd| , ^3| ). The discovery of two configurational tempera- 
tures may suggest the existence of more than one "dy- 
namical temperature" in the off-equilibrium fluctuation- 
dissipation relation approach of [ p~3|JTg| ] . 

We show now how the ensemble averages with the dis- 
tribution (|^) can be analytically calculated for a simple 
toy model of a mixture of two particles. The model can 
be visualized as a cell of the above hard sphere pack, 
made of one large and two small grains seating on the 
vertices of a square lattice, with three layers each of 2 
sites. Interestingly, such a three levels system is able to 
reproduce the general properties found before. In partic- 
ular, the resulting diagram for segregation is plotted in 
Fig.§ in the plane T\ = (mi5a /3i) _1 , T 2 = (m 2 ga (3 2 y 1 - 
It shows that for a given value of T\ by decreasing T 2 the 
system undergoes a crossover, at T 2 *(T\), from a region 
where type-2 particles are on average above (BNE) to 
a region where they are below (RBNE) type-1 particles. 
The structure of the diagram of Fig]4| for small values of 
T\ and T 2 , can be easily understood by looking at eq.(||): 
when the configurational temperature of small particles, 
Ti, is much smaller than the other, T 2 , the potential en- 
ergy of small grains is much smaller than that associated 
to the large ones, therefore large grains segregate on top 
(BNE). The opposite happens when T 2 is close to zero 
and much smaller than T\. In particular, in this toy 
model the RBNE region cannot be met if T 2 is too large. 

The crossover observed in vertical segregation from 
BNE to RBNE corresponds to the drift of the minimum 
of the free energy, F = — InZ, from negative to positive 
values of Ah. In this scenario, the pack undergoing a taps 
dynamics with "amplitudes" Tr (and given "duration", 
t ) at stationarity, as for the data shown in FigJ|, can be 
thus thought to be following a given path, T 2 — T 2 (Ti) 
(where T\ = Ti(Xr, "?o)), in the (T\,T 2 ) plane, crossing 
from the BNE to the RBNE region. 

In conclusion, we have shown that in the present hard 
sphere binary mixture model Edwards' approach to de- 
scribe segregation holds if two thermodynamical parame- 
ters are introduced. In the case considered here these two 
thermodynamical parameters play the role of configura- 
tional temperatures related to the two types of particles. 
By applying all the standard techniques of statistical me- 
chanics of the canonical distribution one can thus predict 
the complex segregation patterns as found above. In par- 
ticular we have shown on a hard sphere mixture on a lat- 
tice that the segregation phenomena in granular media 
can be studied applyingstandard statistical mechanics 
to a new Hamiltonian (|5|), which contains the original 
Hamiltonian of the grains plus a term which constrains 
the system to be in a inherent state. 

In general for a more complex system one might ex- 
pect more constraints to be imposed, leading to more 
than two thermodynamical parameters. In practice, the 
criteria to determine a priori the required parameters can 



be not easily accessible. However more recently we have 
extended data for very low energy [|o| and found that 
only one thermodynamical parameter is necessary to de- 
scribe the stationary state. This seems to be a general 
feature pl| . If this is the Statistical Mechanics 

approach with only one thermodynamical variable may 
be feasible for low energy. 
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FIG. 1. Main frame The relaxation of large grains av- 
erage height, /12, as a function of the number of taps, t, in 
taps sequences with the shown shaking amplitudes, Tr (tap 
duration, to = 10). The superimposed curves are stretched 
exponential fits. Inset A schematic picture of the 3D lattice 
hard sphere mixture considered here. 



FIG. 2. Main frame The difference of the average heights 
of small and large grains, Ah = hi — h-z, measured at sta- 
tionarity is plotted as a function of tap amplitude, Tr. The 
three sets of points correspond to the shown tap durations, 
to. At high Tr larger grains are found above the smaller, i.e, 
Ah < 0, as in the Brazil nut effect (BNE). Below a Tr(r ) 
the opposite is found (Reverse Brazil nut effect, RBNE). In- 
set The Ah data of the main frame are plotted as a function 
of the corresponding average energy, e. The three sets of data 
do not collapse on a single master function, showing that a 
single macroscopic observable, such as e, doesn't characterize 
unambiguously the system status. 
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FIG. 3. Main frame The average density of small grains 
on the box bottom layer, p\, measured at stationarity as a 
function of the height of small particles, hi . Data correspond- 
ing to different Tr and to approximately scale on a single 
master function. The empty circles are the corresponding 
values obtained by ensemble average with the two temper- 
atures Gibbs measure proposed in the text. Upper inset 
The average density of large grains on the box bottom layer, 
p2, obtained for different Tr and to, scale almost on a single 
master function when plotted as a function of the large grains 
height, fi2- Lower inset The average number of contacts be- 
tween large grains per particle, N c , obtained for different Tr 
and ro, scale on a single master function when plotted as a 
function of the system energy, e. 
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FIG. 4. The analytically calculated diagram for segrega- 
tion in the plane Ti = (migao/Ji) -1 , T2 = (m^gao/fe) -1 in 
the three levels model described in the text (/3i is the in- 
verse configurational temperature of the grains of type i) . The 
crossover from Brazil nut effect (BNE) to reverse Brazil nut 
effect (RBNE) corresponds to the drift of the minimum of the 
free energy in the inherent states space, F, from negative to 
positive values of the vertical segregation parameter, Ah. 



